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ABSTRACT 


The problem of determining the stress distribution in a 
plate, which contains straight radial cracks originating from its 
outer boundary or from the boundary of a circular hole inside and 
which IS subjected to a prescribed loading is one of great 
importance. For instance, plane stress analysis of radial cracks 
in an annular disc subjected to centrifugal and/br thermal 
loading is important to extend the use of Fracture Mechanics in 
designing practically important class of structures such as 
bladed disc assemblies, flywheels, clutch plates etc 

One of the most important Fracture Mechanics parameter 
to be determined in such cases is the Stress Intensity Factor 
(S I.F) The present Thesis deals with the determination of 
S I F. for different crack lengths and number of cracks emanating 
from inner/outer surfaces of annular discs of different b/a 
ratios subjected to centrifugal loading by Finite Element Method 
(FEM) The program developed is also capable of finding S.I.F 
for discs subjected to thermal loading. Further, the Temperature 
distribution in a Clutch plate of an Automobile has been found 
using FEM and Finite Difference Methods (FDM) . The concept of 
cyclic symmetry is used to reduce the problem size and Potter’s 
algorithm has been used for the solution Attempts have been 
made also to evaluate S.I.F for rotating annular disc 
experimentally using the Photo-elastic coatings and Strain-gauge 
techniques. 
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CHAPTER 1 


INTRODUCTION 

1 . 1 Introduction: 

A majority of large scale failures an diverse 

structures such as storage tanks, pressure vessels, pipe lines, 
bridges, ships, turbines, generator rotors, flywheels, rocket 
motors, aircraft etc. , have been invariably traced to the 
presence of cracks originating from a site of stress 
concentration These cracks develop and grow to critical sizes 
that bring about complete fracture at loads that were by no means 
anywhere near the critical load of the structure in the uncracked 
condition Fracture Mechanics is primarily concerned with the 
strength of cracked structures or components of a machine It 
provides a methodology by which one can characterise the strength 
of a structure in the presence of cracks 

1.2 The Radial Crack Problem and the Present Work: 

The problem of determining the stress distribution in a 
plate, which contains straight radial cracks originating from its 
outer boundary or from the boundary of a circular hole inside and 
which IS subjected to a prescribed loading is one of great 
importance. For instance, plane stress analysis of radial cracks 
in an annular disc subjected to centrifugal and/or thermal 
loading is important to extend the use of Fracture Mechanics in 
designing practically important class of structures such as 
bladed disc assemblies, flywheels, clutch plates etc. Similarly, 
the plane strain problem is relevant to the study of the action 
of Internal pressure in the interior of hollow cylinders of large 
wall thicknesses with radial cracks originating from the inner 
surface . 


One of the most important Fracture Mechanics parameter 
to be determined in such cases is the Stress Intensity Factor 
(S.I.F.). The problem of cracks emanating from inner boundary of 
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an annular disc with internal pressure has been previously 
studied for large number of cracks [1-7] The present Thesis 
deals with the determination of S.I F for different crack 
lengths and number of cracks emanating from inner /outer surfaces 
of annular discs of different b/a ratios, subjected to 
centrifugal loading by Finite Element Method (FEM). The program 
developed is also capable of finding S.I F for discs subjected 
to thermal loading Further, the Temperature distribution in a 
clutch plate has been found using FEM and Finite Difference 
Methods (FDM) . The concept of cyclic symmetry is used to reduce 
the problem size and Potter’s algorithm has been used for the 
solution. Attempts have been made also to evaluate S I F for 
rotating annular disc experimentally using the Photo-elastic 
coatings and Strain-gauge techniques 

1.3 Thesis Layout: 

Chapter 2 discusses the Finite Element Formulation of 
the problem of Plane Elasticity It discusses the suitable 
choice of the element to model the inverse square root 
singularity required for elastic analysis of crack It also 
describes the method to calculate S.I.F. The Post-Processor 
developed, to present the output data pictorial ly have also been 
discussed S.E.N. specimen subjected to tensile loading has been 
modelled using the crack-tip element chosen and its different 
contours obtained through the Post-Processor have been displayed. 

Chapter 3 explains the Potter’s solution algorithm 
to solve a problem which is cyclically symmetric. 

Chapter 4 discusses the evaluation of S.I.F. in the 
case of a rotating annular disc with radial cracks S I F.’s 
calculated numerically for different cases have been presented 
graphically, ( 03 L“'^ 2 ) » »^r contours developed 
through the Post-Processor have been displayed for one set of 
inner and outer cracks 
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Chapter 5 discusses the evaluation of S.I.F in the 
case of an annular disc with radial cracks subjected to thermal 
loading It also discusses the application of FEM and FDM 
techniques to find the Temperature distribution in an Automobile 
Clutch plate A sample problem has been solved and the results 
are presented graphically 

Chapter 6 discusses the experimental methods using 
Strain Gauges and Photo-elastic Coatings to determine the 
SIF’s of aSEN. specimen and a Rotating disc with radial 
cracks Isochromatics in such cases, observed through the 
Polariscope are also displayed. 
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CHAPTER 2 

FINITE ELEMENT FORMDLATION OF PLANE ELASTICITY 


2 . 1 Introduction 

In this chapter, the general finite element formulation 
for elastic bodies with cracks in a state of plane stress/strain 
subjected to body force, surface force and thermal loading is 
discussed Formulation of the eight noded isoparametric crack-tip 
element is discussed in detail. The value of Stress Intensity 
Factor (S.I F ) is calculated for a Single Edge Notched (S.E N. ) 
specimen by F.E.M. using the crack-tip element. The results 
compare well with analytical solution 

2.2 Problem of Plane Elasticity 

Consider a linear elastic solid of uniform cross- 
section and dimension ’h’ in the z-direction This 3-dimensional 
problem could be simplified to a 2-dimensional problem under two 
different conditions, one, when the thickness ’h’ is very small 
compared with lateral dimensions (say, in x,y directions) and the 
other when the thickness ’h’ is very large These are called 
plane stress and plane strain problems respectively when 
subjected to further conditions as described below ■ 

2.2.1 Plane Stress 

When a thin disc is subjected to support conditions and 
external forces on the lateral boundary which are in the plane of 
the disc and which do not vary across the thickness as shown in 
Fig 2.1, the body is said to be in a state of plane stress 

For plane stress problems, the above assumptions imply 

I 

that all the associated stresses with the z-direction arel 
i V small. i.e.. ' 
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*^22 - "^xz ~ "^yz ~ ® (2 1 ) 

and the remaining stresses are independent of z 
2.2.2 Plane Strain 

When a cylindrical/prismatic body is subjected to 
support conditions and external forces which have no components 
along the axis of the cylinder /prism (in z direction) and the 
remaining components do not change along the axis, as shown in 
Fig 2.2, the body is said to be in a state of plane strain 

For plane strain problems, all strain components 
associated with the z-direction are zero, i.e, 

^22 ~ ^X2 ~ ^yz ~ ® .(22) 

and the remaining components are independent of z. Then, it is 
enough to analyse only a slice of the body 

Although Cg is non-zero in plane stress problems and Og 
is present in plane sti'ain problems, these components can be 
eliminated from the governing equations. 

2.3 Governing Equations of Plane Elasticity Probl^s 

2.3.1 Strain-Displacement relation 

When the deformation is small, the total strain in a 
body in case of plane elasticity problems can be expressed in 
terms of the following strain-displacement relations. 

= 6u/6x, €y = ov/6y, x^y “ 6u/oy + 6v/6x, ..(2.3) 

where u & v are displacement components in x & y directions 
respectively . 



Dsstnbuted force p^, 


Support conditions ^ 

U = L = U 


& 

I =0 I 0 


FIG 2.1 Problem of Plane Stress in Plane Elasticity 



> 



FIG 2 2 Probleni of Plane Strain in Plane Elasticity 


-JU- 
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In the Matrix form, we have 

{€} - [L] {u} (2 3a) 

where the strain vector {€}, the displacement vector {u} and the 
operator matrix [L] are given by 


{€} = - 

1 ^y 

- . CL] = 

O/Ojj 

0 

0 

6/6y 

. {u} = ' 

u 


1 ■’^xy 

L J 

1 

i 

1 

6/by 

6/6x 

1 

V 


(2 3b) 

The total strain at a point in a body is the sura of two 
parts mechanical strain due to stress and thermal strain The 
strain due to temperature change even in the absence of stress is 
called thermal strain. For an isotropic material, symmetry 
arguments show that the thermal strain must be a pure expansion 
or contraction with no shear-strain componenets referred to any 
set of axes. For moderate temperature change, the thermal 
strains due to a change in temperature from Tq to T can be 
expressed as 

= a(T-To) 

= a(T-To) 

= a . .<2.4) 

where a is the coefficient of linear expansion. 

2.3.2 Stress - Strain Relations 


Simplifying the thi’ee- dimensional equations relating 
stress to the mechanical strain for a linear isotropic material 
in the case of plane stress problem (Cg =0), we get 
E 

^x - 2 “ ^xo) ~ ^yo^^ 

1 - V 


o 


y 


E 

2 " ^xo) f^y ~ ^yo^3 

1 - V 


E 


'xy 


2(1 + V) 


^y 


. ..(2.5) 


where ’E’ is the Young’s modulus and ’v’ is the Poisson’s ratio 
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In the Matrix form, these equations can be written as 

{a} = m ({€} - {€> 0 ) ..(2.5a) 

where , 


{0} = 


Ox 


'xy 


{€o> 


IS the stress vector, 


a(T-To) 

cc(T-Tc,) 

0 


is the thermal strain vector 
E 


[D] 


1 - V* 


1 

V 

0 


V 0 
1 0 
0 l-v/2 


(2.5b) 


(2.5c) 


.(2 5d) 


is the elasticity matrix for plane stress 

Similarly, we get the stress-strain relations for a plane strain 
problem (€ 2 = 0 ). These equations can also be put in the Matrix 
form as 


loj = m ({€> - {€o>) 
where , 

E j 1-v 

£D3 = I V 

0 


. . . ( 2 . 6 ) 


{€ 0 } 


(H-v)(l-2v) 
is the elasticity matrix and 


C((T-To) 
(1+v) -I a(T-To) 
0 


V 0 

1-v 0 

0 1-2V/2 


. . . (2.6a) 


. . .{2.6b) 

is the thermal strain for Plane Strain problems. 


2.3.3 Equillbrlun Equation 

Force equilibrium of an infinitesimal particle gives rise to 
the following equilibrium equations 
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OCTx 

+ 

6x 

^'^xy 
+ 

6x 

where 

respectively. 


6y 

OOy 

6y 


+ Fx 
+ Fy 


0 

0 


& Fy are the body forces 


along the x & y 


. (2.7a) 

. .(2.7b) 
directions 


2.3.4 Boundary Conditions 

Formulation of a solid mechanics problem is not 

complete without the boundary conditions The typical boundary 
conditions are as follows First, the total boundary F is 
divided into two parts, i e 

r = Fu U Ft ...( 28 ) 

i) Displacement boundary conditions 

Displacement components are specified on F^ 

*■ 

u = u 

V = v"^ . . (2.8a) 

ii) Stress boundary conditions 

Stress components are specified on F-t 

Ty-ly* ...(2.8b) 

where the surface traction components (Tx.Ty) are related 
to the stress tensor components (cyx.cry, xy) relation 
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*^x "^xy 
'^xy °y 

... (2 8c) 

where and ny are the components of the unit outward normal 

Thus, the formulation of a typical plane elasticity 
problem consists of 

i) Strain-Displacement relations (Eq.(2.3a)) 

11 ) Stress-Strain relations (Eq {2 5a)) 
iii) Equilibrium equations (Eqs (2.7a-b)) 
and iv) Boundary conditions (Eqs (2 8a-b)) 

Even though these equations are linear, they are quite 
difficult to tackle analytically, especially when the geometry of 
the domain is complicated. The alternate method to solve a 
Solid-Mechanics problem is to minimise the total potential 

2 . 4 Total potential energy 



The total potential is given by 

•i fi 

UdV - {F}"^ {u} dV - {T}"^ {u} dA ..(2.9) 

V J V Ja 



The most widely used method for minimization of the 
total potential is the Finite Element method as described below. 
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2.5 8 noded isoparamet;ric element 

The unique description of the displacement within each 
element in terms of nodal values at boundary points or internal 
points of the element is the basic step in any displacement 
finite element formulation and can be expressed as 

' {u} = [N3 {u}*^ .(2 11) 

where {u) is the displacement vector, 

[N] IS the shape function matrix, 

{u} is the nodal displacement vector of the element. 

If the order of the shape functions chosen to describe 
the element geometry are identical to those used to prescribe the 
independent variable, then the element is termed isoparametric, 
A typical isoparametric element is shown in Fig. (2 3 ). The 
functional property of the interpolation function is that 

its value in the natural co-ordinate system is unity at node ’i’ 
and zero at all other nodes For a two dimensional 8-noded 
element, a quadratic variation of the displacement and geometry 
are assumed for which interpolation functions to N3 are given 
in Appendix A Thus, it has 16 displacement degress of freedom, 
two degrees of freedom at each node. Cartesian co-ordinates x , 
y and 2natural co-ordinates r^.s^ are related as follows . 

X = N^xi + N2X2 + . + Ngxg .. ( 2 . 11 a) 

y = N^yi + N2y2 Ngyg ... (2.11b) 

Similarly, the displacements expressed in nodal values are : 


u - Njui + N2U2 + . + Nsub 
V = Nj^V]^ + N 2 V 2 + + N0V8 


. . (2,11c) 
.. (2. lid) 
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s A 

I 



X 


F\g.2^ Isoparametric transformation fora quadratic element 
from global to normalised co-ordinates. 





X 

Fig E *4 2-D triangular element with mid sides rK>des at quarter 

point from rectangular isoparametric element. 
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Therefore, the sfiape function matrix [N] and the nodal 
displacement vector {u}^ can be written as 


(2 12a) 


Ni 0 N2 
0 Ni 0 


. .Ns 0 
0 Ni 


{2.12b) 


When the discretisation (2.11) is substituted in the 
Eq (2 9) for total potential, we get 

Total potential 


1=2 (5g{u}^ CK]^{u}® - {u}® {to - {fb> ) - 2 ({fs> {ur) 

e=l b=l 

. {2.13a) 

where n© is no. of area elements & n^, is no. of 


boundary elements & 

[K]^ = element stiffness matrix 


[B3 tD][B3 t dxdy 


. .{2 13b) 


{f^}® = Element Body force vector 


[N] {F} t dxdy 


. . . {2.13c) 


{f-t)^ = Element Thermal force vector 


t dxdy 


. . (2 13d) 
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&. = Element surface force vector 

% 

= [N]'^{T} t db . (2 13e) 

i 

The above expression for I can be wiitten in terms of 
the Global Matrices as 

' 1 = ^ {U}'^[K3{U} - {Fb}'^{U} - - {Fs}'^{U} 

- .. (2.14) 

where {U} is the Global displacement vector and the 
Global Matrices [K] , {Fbl. {F-t}> {Fs) contain the contributions 
from all the elements Minimization of the expression (2.14) 
with respect to {U} gives, 

[K3{U} = {F} .(2 15) 

where {F} = {F^} + {Ft) + {Fg} (2 15a) 


2.6 Crack-tip Element Formulation 

2-6.1 Elastic Singularity at the Crack-tip 

Linear elastic fracture mechanics is based on the 
concept that the stress tensor (Oj^j) near the crack tip (taken as 
origin of the polar co-ordinate system) is given by the 
Westergaard solution in the form given by Irwin (Broek [13). 

Oij = K (27tr)"^ fij(e) ..(2.16) 

where, K is the Stress intensity factor and designated as Kj, 
Kjj, and Kjjj for the fracture modes I, II, & III respectively 
fj[j(8) represents the distribution of stress and depends on the 
mode of fracture. 

The above expression shows that at the crack-rip, 
singular stress field exists and it is popularly known as inverse 
square root singularity. 



Aliy finite eleiiient piocedure siiould take into account 
this fact while modelling the crack tip Earlier investigators 
used CST (constant strain triangle) element of very small size of 
the order of 1/100 (1 is tlxe crack length) near the crack- tip to 
model the strain gradient appropriately However, in the last 
few decades many new elements have been proposed by different 
investigators which take care of the inverse sequare root 
singularity and one can afford to use larger size elements near 
the crack tip Elements based on classical solution displacement 
fields, polynomial displacement fields and Isoparametric concepts 
as well as hybrid formulation are reported in the literature 

Afliong the above elements quarter-point singular 
elements are very popular essentially due to the availability of 
the quadratic isoparametric serendipity element in almost all 
general purpose programs 

2.6.2 Historical Overview of Quarter Point Elements 

Barsoum [8] and Henshell and Shaw [9] have 
independently demonstrated that the inverse square root 
singularity characteristic of linear elastic fracture mechanics 
can be obtained in the two-dimensional 8-noded isoparametric 
element (Q8) when the midside nodes near the crack tip are placed 
at the quarter-point. This concept was subsequently extended to 
plate bending and shell fractures by Barsoum [10]. 

Barsoum [11] then showed that the triangular element 
formed by collapsing one side of the Q8 led to far better results 
than the rectangular element. Freese and Tracey [12] have shown 
that the natural isoparametric triangle (NIT) and the collapsed 
quadrilateral perform equally well If the side opposite to the 
crack tip is curved, then there is a substantial deterioration in 
the SIF calculation from the collapsed quadrilateral, and no 
change in the NIT. Such a discrepancy is caused by the fact that 
a collapsed 8-node quadrilateral does not actually degenerate 
into a NIT. 



15 


The extension of the quadratic isoparametric quarter- 
point element to cubic isoparametric was proposed by Pu et al 
[13] Yamada et al [14] extended the concept of the 6-noded 
isoparametric element to the variable numbering element 

Hibbitt [15] proved that the singular rectangular 
element has a singular stiffness while the triangular one does 
not, and that in the collapsed quadrilateral element the 
singularity prevails along the two sides only, while it is 
omnipresent inside the triangular element He attributed to this 
difference the better results achieved by the triangular element 
as reported by Barsoum [11] He also demonstrated that a variety 
of stress singularities (1/n) can be achieved by isoparametric 
elements with a polynomial approximation of order n In a recent 
paper Ying [16] has shown that in his investigation Hibbitt [15] 
erroneously concluded that the strain energy of a rectangular 
quarter-point element is singular, but that, as previously known, 
the singularity is along the edges and diagonal only. Also 
discussed by Ying [16] is the error associated with the location 
of the ’quarter point’ of a singular element This was further 
elaborated by Barsoum’ s [17] discussion of the previous work, 
where he has shown that the error in estimating the SIF is one 
order of magnitude smaller than the discretization error in 
locating the ’quarter point’ Hence, the exact location of the 
roidside node is not crucial as long as the discretization error 
is small. Also, by having multiple independent nodes ( of a 
collapsed Q-8 element), Barsoum [18] has shown that small-scale 
yielding (characterized by 1/r singularity) could be modelled 

Lynn and Ingraffea [19] generalized the concept of the 
quarter-point singular element, and have shown that by varying 
the placement of the side node, between quarter and mid-point, 
one can control the point where singularity is to occur (between 
the corner node and infinity, respectively). This led to the 
introduction of the transition element which, when inserted 
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around the singular elements, resulted in improved SIF 
calculations as 1/a ratio is decreased. Again, Hussain et. al 
[203 extended the concept of quadratic transition element to 
cubic transition elements 

Lynn and Ingraffea [19] investigated both the effect of 
the 1/a (element length over total crack length) and aspect ratio 
of the singular elements. An optimum ratio of 1/a = 0.25 was 

reported, while it was found that the aspect ratio effect was 
relatively unimportant Another detailed investigation of the 
optimum quarter-point element size was carried out by Ingraaffea 
and Manu, [21] where 1/a was varied from 0.2 to 0 03 in both two 
and three-dimensional analyses The errors in the two- 
dimensional results varied from 8 percent for 1/a =0.2 to -1 
percent for 1/a = 0 03 

An assessment of the quarter-point elements (in 
pressure vessel fracture analysis) was offered by Barsoum [22] 
for both two and three dimensional analysis Three sources of 
modelling errors, associated with the quarter-point elements, 
were discussed: (a) those associated with the type of quarter- 
point element, triangular or rectangular, (b) those associated 
with the configuration of the element boundary, straight or 
curved; and (c) those associated with the location of the 
’quarter point’ node. Barsoum [22] indicates that (a) triangular 
elements are to be preferred over the quadrilateral ones (whether 
collapsed or not); (b) the sides of the element should not be 
curved; and (c) a perturbation of the ’quarter point’ node by e, 
leads to an error in calculating the stress intensity factor of 
ge, where g is the ratio of crack tip element to crack length 

Harrop [23] qualitatively discussed the optimum size of 
quarter-point crack tip elements He indicates that the singular 
element can represent the stress singularity and a constant 
finite stress term only Thus, any singular element which 
is ’too large’ cannot represent nonlinear (and non-singular) 
stress variation in a structure On the other hand, by using a 
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singular element which is ’too small’, the error in representing 
the finite stress term decreases, but the region of the mesh 
representing the stress singularity also decreases This source 
of error would still hold even if transition elements were used. 
He thus pointed out that some crack tip element size has to be 
optimum, and concluded that ’ . it is clearly impossible to 

recommend a particular crack tip element size suitable for all 
situations. ’ 

Saouma and David Schwemmer [24] conducted a parametric 
numerical evaluation of the ’quarter-point’ singular elements for 
272 test cases The modified parameters were the order of the 
integration, aspect ratio, number of elements surrounding the 
crack-tip, use of transition elements, the singular element 
length over the total crack length, and the symmetry of the mesh 
around the crack tip Based on their analysis the following 
recommendations were made for the ’best usage’ of the quarter- 
point singular element It is claimed that the SIF calculations 
are within 10 % of the anal 5 ^ical result 

1 Use (2x2) reduced integration scheme. 

2- Use at least four (in pure mode I problems) singular 

elements around a crack tip 

3. Have the internal angles of all the singular elements, 
around the crack tip approximately equal to 45 degress. 

4 Unless an excessive small 1/a ratio is used little 
improvement is achieved by using transition elements. 

5 For problems, with uniform non-singular stress distribution, 
little Improvement is achieved by using small 1/a. 

6 For problems where a non-singular stress gradient is 
expected, 1/a should be less than 0.5. 

These recommendations were taken as guidelines for using the 
Q-P elements in the present work. 
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2.6.3 Quarter Point Element Implementation Consideration 

Singular triangular element is obtained by collapsing 
the nodes 1,4 and 8 and moving the nodes 7 and 5 to the quarter 
point as shown in Fig (2 4). Thus, globally the nodes 1,4 and 8 
are assigned the same co-ordinates and the same node number 
Barsoum [18] has shown that by assigning different node numbers 
for nodes 1,4 and 8 globally, one can simulate 1/r singularity 
which IS prevalent in elastic-perf ectly plastic problem 

2.6.4 Stress Intensity Factor Calculation 

One of the most important param eter to be determined 
in fracture mechanics analysis is the SIF It is a function of 
the loading on the cracked configuration and of the size and 
shape of the crack and other geometric boundaries It has the 
dimension of stress * 4( length). Crack tip opening displacement 
(CTOD) method has been used to evaluate SIF In CTOD method, 
displacement Uq at the quarter point of the singular element on 
the face of the crack is noted and the following relation [6] is 
used to calculate the SIF. 

E U0 1/2 

Ki = (27r/r) . .(2.17) 

(1+v) (1+K) 

where , 

r IS the distance from the crack tip, 
and K = (3-v)/(l+v) for plane stress and 
(3-4v) for plane strain 

2.7 Post Processor [253 

Postprocessor displays the output data in a meaningful 
way The deformed structure can be displayed, either 
superimposed on its undeformed state or contours of stress, 
displacement can be plotted. 
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The procedure for plotting the contours is simple 
Each element is davided into a number of triangles (in the 
present case, six) such that vertices of the triangles are the 
nodes. The element is so divided that the triangles do not 
overlap each other The program assumes linear variation along 
the sides of the triangle Depending upon the value for which 
contour is plotted, the point is located on each sides of the 
triangle, if present. These points are joined together to form 
the contour in that triangle The program also has the fill 
option. In this, the zone between two successive contours is 
filled with a suitable colour The program developed can use 256 
different colours and the range of each coloured contour is also 
displayed. 


The plotting of displacement contours is straight 
forward since in the displacement formulation of FEM, the results 
obtained directly give these values at the nodes. For plotting 
contours of stress, one has to evaluate these, based on 
displacement values In some cases transformation to global 
cartesian will be required 


As stresses are not accurate on the element boundary 
[263, they have been calculated at the gauss points of each 
elements. These stresses are extrapolated to the corner nodes by 
the following relation, 

Me = CT] [o3g 

where [ol^ is stress at corner nodes, 
lalg is stress at gauss points, 

and 


[T] 


l+(^3/2) -0.5 l-{*f3/2) -0 5 

-0,5 l + (43/2) -0.5 l-('f3/2) 


l-(-r3/2) -0.5 l+('r3/2) -0.5 

-0.5 l-(^3/2) -0.5 l+(^3/2) 


is the transformation matrix 

The values of stress components at midside nodes are averages 
of the values at the corner nodes adjacent to the midside nodes 
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i.e , 

as = {ai ■+ a2)/2 

0Q = {02 + cf3)/2 

<^7 = {02 + a 4)/2 

as = {a4 + ai)/2 

The nodes which lie on more than one element will have 
different values of stress corresponding to each element The 
value averages over all the elements to which the node belongs is 
considered as the stress at that node. 

2.8 EVALUATION OF S.I.F. IN AN S.K.N. SPECIMEN 

In order to test the 8-noded isoparametric 
quadrilateral crack-tip element, the problem of Single Edge 
Notched (S E.N. ) specimen subjected to tensile loading is taken 
and its S.I F. (Kj) is evaluated The mesh chosen for the problem 
is shown in Fig (2.5) 

The stiffness matrix discussed in Sec. (2 5) is referred 
to cartesian co-ordinates Due to symmetry of the problem, at 
*'■*-21 & 267-272, Y- displacement is suppressed and at nodes 
253-268, X-displacement is suppressed. To impose the Force 
boundary condition, the product of the applied stress (a) and 
element length on which the stress is applied, is distributed in 
the ratio 1/6 : 1/6 • 2/3 for two corner nodes and mid-side 

node respectively 

The assembled stiffness matrix in banded form is solved 
by the Gaussian Elimination method. C T 0 D method has been used 
to calculate the S.I.F. 

Analytically, the value of S.I.F. for S E.N. specimen 
is given as, 

Ki = ( 1.1215 - 0.23 (a/w) ) o -fxa 

where a/w is the ratio of the crack length to width of the 
specimen 

The results from both the methods are compared in 
Tab. (2 1), Various contours such as {ai-a 2 ) . (^ 1 + 02 ) .Ux.cry , 
Ux,Uy and deformed mesh, developed through the Post-processor are 
shown in Figs (2.6-2 11). These figures validate the 
correctness of the crack-tip element used 
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FIG.2-5 MESH FOR S.E.N. SPECIMEN 
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2.9 Closure 

In this chapter, Finite Element Formulation of Plane 
elasticity problem subjected to Body force, Surface force and 
Thermal loading has been discussed. Also modelling a crack, with 
8-noded isoparametric quadrilateral element is discussed It is 
seen that a large number of elements are required to model a 
crack-tip In the case of a radial crack problem, if one solves 
a problem as a whole, the number of elements will phenomenally 
increase as the number of cracks increases. However, if one 
takes into account the symmetry of the problem, the amount of 
computation can be reduced. The exact requirement of incore 
memory and CPU time largely depends upon the solution procedure 
Next chapter discusses Potter’s solution algorithm which uses 
minimum amount of incore memory for solving a cyclically 
symmetric problem. 
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CHAPTER 3 

SOLUTION ALGORITHM FOR SOLVING A CYCLICALLY SYMMETRIC PROBLEM 

3 . 1 Introduction : 

Annular disk without crack can be easily analysed as an 
axissnjimetric problem But as soon as cracks are formed the 
structure looses its axisymmetric nature. Modelling of a crack 
involves the use of a large number of small elements near the 
crack tip. As the number of cracks increases, analysing the 
problem as a whole becomes prohibitive and one would require 
large computer resources to solve the problem However, if the 
cracks are periodic and of equal length, then this periodicity 
can be utilised in solving the problem The assumption of 
periodic cracks of equal lengths has been used by all earlier 
investigators [3-6] and this assumption has been experimentally 
verified by Shukla [7] for the case of internal pressure loading 
However, extending the same logic one can assume that cracks 
occur with certain periodicity in rotating discs too. In the case 
of Clutch plates in an Automobile, failure has been observed with 
three radial cracks emanating from inner boundary. 

3.2 Solution Algorithms for a Cyclically Symmetric Problem - 

A problem is cyclically symmetric when geometry of the 
structure, material properties, loading and boundary conditions 
are also cyclically synxnetric. Since, the cracks occur with 
certain periodicity, the problem of an annular disc with radial 
cracks can be idealised as a cyclically symmetric problem 

In the case of static analysis, the displacement vector 
of each sector are identical and it is enough if one analysed 
only one repeating sector of the whole structure. In order to 
impose the condition that the structure is cyclically symmetric, 
it is desirable that the globixL stiffness ma, trix is referred to 
the global polar co-ordinates. 



3.3 Transformation from Global Cartesian to Global Polar : 


For a cyclically synuaetric problem, to apply the condition 
that the displacement vector of the first and last partition is 
same, it becomes necessary to convert the problem form global 
cartesian co-ordinate system to global polar co-ordinate system. 


For any node point shown in Fig. (3.1) the ti-ansf ormation is 
defined as 


Xl 

- = [Ti3 ■ 

^1 

Wliere Tj, = • 

COS^j^ 


Yl 


©i 


sin^j^ 

COS(6j^ 


— (3 1) 



I'l 

T 

r 

X 

or ' 

01 

- = CTi3 - 

y 


and = tan ^ (Yi/xj^) 


Therefore, for an element, co-ordinate transformation is defined 
as, 

{U} r [T3{U}p (3.2) 

{F} = CT3{F}p 

where the subscript ’p’ refers to the polar co-ordinates and the 
Matrix [T3 has on its diagonal the submatrices [Ti3 depending on 
the co-ordinates of each node Pre-mult iplying the expression 
(3 15) by [T3'^ and using the Eq (3.2), we get 


i;K3p{0}p = {F>p ...(3.3) 

where 

[K3 = CT3^[K3 CT3 

' • ^ J P 

Eq.(3.3) can be solved by many techniques. One such 
efficient method viz., Potter’s method which has been used to 
solve the equation is explained below. 
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3.4 Potter’s Method £273' 

This IS a method of substructuring The basic approach is 
to split the problem into a number of small problems The 
concept of partitioning helps in understanding the method. 

The stiffness equations for a typical substructure with 
extra partition Fig. (3 2) is, 

[K]{U} = {F} (3.8) 

The linear problem such as shown in Fig (3 3) will have 
global stiffness matrix in banded form This problem can be 
divided into four substructures I, II, III & IV are the 
dividing lines of the substructures & are known as partitions. 
Each partition is identified with the number of nodes it 
contains 

Partition I contains node no. from 1 to 4 partition II, 
contains node no from 5 to 8 so on. 

Partition is nothing but a convenient way of representing 

the global stiffness matrix. Effect of partitioning on global 

stiffness matrix is shown in Fig. (3.4). Matrix and 

are the submatrices and their meaning is as follows. 

[B^] - It represents the effect of nodes of i partition on 

itself . 

CAi3 - It represent.. the effect of nodes of (i+l)^ partition 

on nodes of i^ partition 

CCi3 - It represent.j^the effect of nodes of (i-1) partition 

on nodes of partition 

Equation (38) can now be represented as. 


gl Ai 


Zl 


Si 

^2 ®2 A 2 


22 


S2 

C 3 B3 A 3 


Z3 



C 4 B 4 

- « 


Z4 


G 4 


(3 9) 
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34 ■' ? 

Where, {Z^} is the displacement vect^ of i^^ partition 
{G^} is the force vector of i^^^ partition 

Potter’s Method does not calculate displacement of nodes 
individually, but instead calculates displacement vector for each 
partition (i.e ) Zj^’s This means that for example when it 
obtains the displacement vector of partition I, displacement of 
node 1 to 4 are obtained simultaneously. 

Potter’s Method is an out of core solution technique At a 
time it requires only submatrices of the i^ partition in the 
core memory of the computer All other matrices are stored on 
secondary storage devices such as taper, disks This makes 
optimum use of the computer in core memory 

3,4,1 Global Stiffness Matrix and Sector Stiffness Matrix for 
A Cyclically Symmetric Structure: 

In a Cyclically symmetric structure, a substructure can be 
identified as the one which repeats itself periodecally. 


Consider a symmetric structure as shown in Fig (3 5) It 
has H substructures and each is divided into n partitions, thus 
the total number of partitions are m (=Nxn) Unlike in Fig 
(3.3), here the last partition nemely the partition and the 
first partition are connected. Hence, by the definition of 
submatrices [A], [B] and [C3, in the present case matrices and 
Ajn exist. Thus the global stiffness matrix takes the form, 



Al 

0 

Cl 

C2 

B2 

A2 


^m 

0 


Cm Bffi 


(3.10) 


Partition 1 & (n+1) constitutes as sector boundaries. Since the 
structure is cyclically symmetric, in the static case the 
displacement vector of each sector is same i.e.. 


{Zn+l} = {Zl). {Zn+2) = W 

{Zn+j} = {Zi}. {Z2n} = tZxi) <3-1“ 

Osing the above condition the overall problem can be reduced to 
only one sector and the stiffness matrix takes the form : 
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Bl 

Al 

Cl 




G2 

B2 

A2 



. (3.10) 

Am 


Cn Bn 




Due to 

finite element formulation of 

the problem. 


ECi3 ^ 

= [An 

3*^ and [C^+i) 

= [All" 

for i = 1,2, 

, n-1 

3.4.2 

Evaluation of Displacement vectors: 



The static equilibrium at each station can be represented as 
follows . 

For the first station, 

(Zn) + [Bi3 {Zi} + [Ai] {Z 2 } = [Gi3 (3.14) 

Similarly for 1 ^ station, 

CCi3 {Zi-i} + [Bi3 {Z^} + [A 3.3 {Zi+il = [GO (3.15) 

For the n^ station, Aj-^ the coupling stiffness matrix with the 
first station of the next sector The above equation can be 
modified as, 

{Zi} = -[Pi3 {Z 2 } + [qi3 + [Qi3 {Zn> (3.16) 

and 

{Zi} = - EPi3 {Zi+i) + [qi3 + CQi3 (Zn) (3.16-A) 

where , 

[Pl3 = [Bi"^3[Ai3 
[qi3 = EBr^3[Gi3 
[Ql3 = [Br^3[Ci3 
CPi3 = (CBi3 - CCi3Ci-i3)'^ [Ai3 
Cqi3 = ([Bi3 - [Ci3[i-i3)"^ ({Gi) - CCi3{qi-i}) 

[Qi3 = -([Bi3 - [Ci3[i-i3)“^ ([Ci3[Qi-i3) i = 2,3, — n (3.18) 

Considering the equilibrium at the n station, 

tCn3 {Zn-l) + [Bn3 {Z^} + [Ai^l iZn+l) = [Gn3 (3.19) 

Because of Cyclic syinmetry {Z^+i) = [Zi3 

Hence , 

CGn3 {Zn-l) + [Bn3 CZn) + CAnD {Zi} = [Gn3 


(3.20) 
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Substituting for {Z^} from Eqn (3.16) gives, 

[Cn3 {Zn-i) + [Bnl ■+ [A^] ([~Pi]{Z2} + + [QiliZ^}) 

= t<5n3 (3.21) 

1 e 

CRl] {Z2} + [Si3 {Z^} + [Ci,3 {Zn-i) - [ui] (3.22) 

where, [Rll = [-Aji3[Pi3 

[Si3 = [Bn3 + [An3 [Qi3 
[ui 3 = [Gn 3 - [An 3 [qi 3 
{Z 3 }, ---{Zn- 2 > gives, 

[Ri3 {2i+i} ■+ [5i3 {Z^} + [C^] (Zn-i) = [u^] (3.23) 

Here, [Ri3 = i:-Ri-i3[Pi3 

tSi 3 = ESi_i 3 + [Ri-i 3 [Qil 

Cui3 = [Wi-i3 - [Ri-i3 [qi3 1 = 2,3 n-2. 

Rearranging the (n-2)nd equation yields. 

([Rn-23 + EBn3{Zn-i} + ([Sn-23{2n} = {'in-2> (3.24) 

Substituting for {Zn-f} from Eqn (3 10) gives, 

[([Rn-23 + [Cn3) ( [Qn-1^ ~ fBn-l3) + CSn-23 3 {Zn) 

= [{Un-2> - ([Kn-23 + [Cn3) (qn-l) (3.25) 

After evaluating {Zjj} a series of backward substitutions in Eqs 
(3 16) & (3.17) will provide all the unknown displacements 
■t2n-l3^> { 21^-2 {Zj} 


3.5 Closure : 

In this chapter. Potter’s solution algorithm has been 
discussed to solve a problem which has cyclic symmetry. In the 
next two chapters, the radial crack problems subjected to 
centrifugal and thermal loading have been extensively analysed 
using the software developed in this chapter. 



38 


CHAPTER 4 

E7ALDATI0N OF S.I.F, FOR ROTATING AHNGLAR DISC WITH RADIAL CRACKS 


4.1 Introduction : 

In this chapter, the results of the analysis performed 
on a rotating annular disc with radial cracks using the software 
discussed in the previous chapter is presented S.I.F s have 
been calculated as a function of crack length, b/a ratio and 
number of cracks. First S.I.F s as function of number of cracks 
(n) are calculated for three crack length ratios with {l/(b-a)} 
as parameters Such graphs are obtained for two different b/a 
ratios Next, S I F s as a function of crack length are 
calculated for three b/a ratios as parameters and such graphs are 
obtained for two different number of cracks. These calculations 
have been done for both inner and outer cracks. 


4.2 Evaluation of Elraaent Body Force Vector : 


In the previous chapter the element body force vector 
has been derived as; 


{fb) 


[N]^ {F} t dxdy ...(4.3) 

The vector {F} can be expressed in teims of nodal values as, 

{F} = 


0 H2 ® 

0 Nl 0 N2 


Nb 0 

,0 Ng 


Fxl 

Fyl 


Fx8 


= [N3 {Fb}' 


(4.2) 
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Thus, 

= 1 i [N]" m t dxdy {Ffe}^ ....(4.3) 

Nuw, the centrifugal force at a point in polar co- 
ordinates IS given by, 


Fr = 


. (4.4) 


where , 


= density of the material, 
w = angulr speed of rotation and 
ri = radius at node i 

Transforming into cartesian co-ordinates, the nodal 

forces are given by, 

- Fj^^cos 

Fy^ = Fr^sin d± 

•where 6]^ = tan ^iy±/xi) .. .(4.5) 


This integral (Eq (4 3)) when evaluated by Gauss-quadrature 
technique yields the consistent loads at each node as shown in 
Fig. (4.1). Looking at the Fig, (4.1) , one is surprised to see 
that the nodal forces are unequal and the corner nodal forces 
are negative However, the mid-side nodal forces are positive and 
in the direction of body force. The apparent ambiguity is due to 
consistent force formulation, A similar effect is explained by 
Robert D.Cook [283 when plane 8-noded quadratic element is 
chosen. Further, he has shown that the nodal forces are equal and 
in the same direction as of body force when the plane 4-noded 
linear element is chosen. 


4.3 Preparation of Input Data : 

As large amount of input data is required to solve the 
problem, it is difficult to prepare the complete input file 
manually. Instead programs have been prepared to aid in 
preparing the data. The procedure adopted is explained behind. 
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To prepare input file for a problem, it is divided into 
3 parts as shown in Fig. (4 2). 

The program basically prepares co-ordinates of each 
node and nodal connectivity of the eight noded isoparametric 
elements for a sector of an uncracked disc The input data for 
the program is number of radial lines, number of divisions in the 
radial lines, angles at which radial lines appear, radial 
distance of the radial divisions and the node number at which 
first node will start. This program is sufficient to generate 
input data for part 1 and part 3 mentioned above 

For part 2 the same program is used initially to 
calculate nodal co-ordinates and nodal connectivity as though it 
has no crack tip. Later these output become the input for 
another program via , ’inn f’/ ’out f’ depending on whether it 
is an inner crack problem or outer crack one, which introduces 
the Quarter Point (Q-P) elements at the crack tip 

Thus, after getting the modified output file from 
’inn.f’/ ’out f’ for part 2, it is added with the output files of 
part 1 and part 3 suitably to generate the complete input file. 

4.4 SuBiierical Results : 

In order to obtain accurate results, close partitions 
are formed near the crack line and it becomes less dense as it 
goes away from the crack line as shown in Fig (4.3) and 
Fig, (4 4) As the number of cracks increases the sector angle 
becomes smaller and smaller. Thus contrary to the total size of 
the problem the CPU time and memory decreases as number of cracks 
increases. This is a tremendous advantage as fracture mechanics 
problem requires large number of elements in the vicinity of the 
crack tip. The S.I.F. is calculated using CTOD method. In all the 
cases, 9 displacement is suppressed on sector boundaries. It is 
justifiable because of symmetry. 





yiG 4 2 Finite Element Mesh showing Partitions 
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N 


FIG 4 3 Finite Element Mesh for a sector with one 

Radial crack originating from Inner boundary 
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Crack tip surrounded 
by Q-P elements 

Crack line 


FIG. 4. 4 Finite Element Mesh for a sector with one 

Kadial crack originating from Outer boundary 
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4.4-1 Periodic Redial cracks originating from inner boundary : 

Firi,i. n. J.F.s as a function of number of cracks 
(n)=:s4,?,12,::0.::4 h 36 are calculated for three crack length 
ratios namely 0. 1 , 0.5 & 0.7 and each for two different 

b/a ration vi/. . <•. h. 10. Next, S.I.F.s as a function of crack 
length raiio.s (3/h~a) from 0.1 to 0.8 in steps of 0.1 are 
calculated for three b/a ratios namely 2,3 & 10 and each for two 
different number of cracks via. 8 & 4. The results are non- 
dimensional Ised by dividing it by =( (3+v)/8)*( f w^(b^-a^) ) . 
These results are shown graphically from Figs. (4. 5-4. 12) . 

It is observed from the graphs that, 

i) .j. I.F. decreases as the number of cracks increases 
for constant crack length and b/a ratio. 

ii) i,i . I . B . increases as the crack length increases for 
the same number of cracks and b/a ratio. 

iii) S.I.F. decreases as b/a increases keeping b 

constant fox* the same crack length ratio and number 
of cracks. 

The first two observations are in line with common- 
sense, but the third observation is less obvious. In changing 
b/a ratios, the inner radius was reduced to increase the b/a 
ratio. Fig. (4. 13) shows the variation of og and Oj- at the 
central radius of an annular disc (without crack) as a function 
of b/a ratio. In this case, the variation of b/a was achieved by 
keeping ’b’ constant and ’a’ was reduced. Fig. (4. 14) shows a 
similar graj>h wherein ’a* was kept constant and ’b’ was increased 
to get various b/a ratios. Since in an uncraoked annular disc, 
00 decreases when ’a’ is reduced, the S.I.F, value also 
decreases. If b/a ratio was changed by increasing ’b’ then, 
S.I.F. would have increased for increasing b/a ratios. This 
shows that the b/a ratio alone is not sufficient to represent the 
results in non-dimensional form. 
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4.4.2 Periodic Radial cracks originating from Outer boundary : 

S I F b are calculated similarly for all the 
parameters ai that uf inner cracks and the results have been 
shown grapiiically from Figs (4.15-4.22). 

It IS observed from the graphs that similar trend is 
present in outer cracks as that of inner cracks (Ref Sec 4 3 1) 
except that the magnitude of S I.F. is little greater for an 
inner crack as compared to an outer crack of the same dimensions. 

4.4.3. Display of output data : 

In order to validate the correctness of the solution 
algoritlim, vari<»us contours such as (01-02). (01+02). (o^). (09), 
(Ur)> (1^9) Deformed mesh are drawn using the Post-processor 
developed These are shown in Figs. ( 4. 23-4. 36) . The symmetry 
of these indicates quail tatively that the solution algorithm is 
correct. Figs. (4.37 & 4.38) show the Mode-I isochromatics 
(contours of 01-02) observed in a Photo-elastic experiment [28] 
for annular discs subjected to internal pressure with cracks 
emanating from inner and outer boundaries. Though in the case of 
rotating discs, the loading pattern is different, only Mode-I 
exists at the crack -tip due to centrifugal forces The contours 
of (01-02) for both inner and outer cracks are qualitatively 
similar to Figs. (4 . 37-4 . 38) . This completely validates the 
correctness of the algorithm used. 
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4.4.4 Plastic zone at the crack-tip : 

In this chapter, only Linear Elastic Fracture Mechanics 
(LEFM) approach is used for evaluating S I F. In actual 
engineering materials, plastic-zone invariably exits at the 
crack-tip and if the zone is small compared to the crack length 
then, LEFM approach is valid Plastid zone is evaluated by using 
the Von-Misces yield criterion as given below; 

2 ^ 2 2 
ai +02 

In the FEM modelling, disc was considered to be of 
Aluminium of the values of E 70 GPa and Oyield - 
plastic zone for various RPM are obtained for both inner and 
outer cracks. From the Figs. (4 39-4 40), it is clear that the 
plastic zones are very small and LEFM approximation is valid 

4.5 Closure : 

It is successfully demonstrated that the concept of 
cyclic symmetry is very helpful in determining the S I.F in the 
case of a rotating disc with radial cracks. In the next chapter, 
evaluation of S.I F in an annular disc subjected to thermal 
loading IS discussed. 
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FIG 4 13 Stresses at the Central radius of a Rotating 
Annular disc for various values of b/a 
( ’b’ increases and ’a’ is constant ) 
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FIG 4.14 Stresses at the Central radius of a Rotating 
Annular disc for various values of b/a 
{ ’a’ decreases and ’b’ is constant) 
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FIG, 4 23 



FIG. 4.24 
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FIG. 4.26 
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FIG. 4.32 



MB4t 


INVESTGN. OF OITER CRfiCK 


III-B 




FIG. 4. 38 foj-02> coratours obtained experimentally in an 

Annul ar disc witli 4 Outer cracks subjected to 
internal pressure [28] 
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CHAPTER 5 


EVALUATION OF S.I.P. OF AN ANNULAR DISC 
SUBJECTED TO THERMAL STRESS 


5.1 Introduction : 

In this chapter, evaluation of S.I.F. in an annular disc 
subjected to thermal loading is discussed. It has been observed 
that the Clutch plate in an Automobile fail due to cracks 
emanating from inner boundary. An attempt has therefore been made 
to find the temperature distribution in a clutch plate enabling 
to calculate the S.I.F. of cracks in such plate. 

5.2 Evaluation of Element Thermal Force Vector : 


In Chapter 2, the element Thermal force vector has been 
derived as: 


(ft) = 


I 

i 


CB3^ tD}£€o} t dxdy 


...(5.1) 


But {€o> « CT) 

Thus, if the Temperature distribution (T) in a disc is 
known, the stresses due to thermal loading can be evaluated by 
FEM. 


5.3 Wume rioal Besulta : 

A sample clutch plate of certain dimensions is taken 
and its Transient Temperature distribution is obtained by 
coidjiaation of FEM and Finite Difference Methods (FSM), The 
solution procedure is attached in Appendix . The Teaperature 
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distribution sc> ..btainrd is shown in Figs. (5. 1-5. 3). ia^-C2) and 
value.s hav«- also been calculated for such plate and these 
contours are shown in Figs .( 5. 4-5. 5) . oi being Tnaxinum at the 
inner boundary possibly explains the cause for cracks emanating 
from the inner boundary in clutch plates. 

If such a clutch plate has 4 cracks from inner bounday, then 
its S.l.F. has been evaluated using the software developed by 
Cyclic Symmetry concept. The results have been obtained 
considering different temperature distributions such as T = r, T 
z r^i T = l/r, T ^ log r, T = log r + 503 and the actual 

distribution as obtained above. The results have been tabulated 
in Tab, (5.1). The (oi-02) contour developed for logarithmic 
distribution is shown in Fig. (5.6). The ( 01 - 02 ) contour compare 
well with the isochromatics corresponding to Mode-I type of 
loading. This qualitatively indicates the correctness of the- 

implementation . 

5.4 Closure : 

It has been shown that the concept of cyclic symmetry 
can be used effectively to determine the S.l.F. in the case of an 
annular disc subjected to Thermal loading. Hext chajpter deals 
with the Experimental Methods to determine the S.I.F.’s. 
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FIG. 5. 2 Transient Temperature distribution in a Clutch 
plate at the Outer Radius as a function of 
Distance in Axial direction 
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FIG. 5 3 Temperature Distribution in a Clutch plate 
at the Surface as a function of Radius at 
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S No 

Distribution 

S I F 
{ MPa-fm ) 

1 

T = r 

0.129 X 10~^ 

2 

T = r® 

0.626 X 10'^ 

3 

T = 1/r 

0.224 X 10'^ 

4 

T = log r 

j 0 536 X 10~^ 

5 

T = log r + 500 

j 0.669 X 10~^ 

6 ■ 

as obtained in 
tlie analysis 

j 0.420 X 10"^ 


51 S.I F’s for a Clutch Plate with 4 Inner Cracks 
Suhoected to Various Thermal Loadings, 


Table 
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CBAPTER 6 


EXPERIMENTAL DETERMINATION OF STRESS INTENSITY FACTORS 


6 . 1 Introduction : 

In the previous chapters, the use of Finite Element 
Method in the determination of S I F for the radial crack 
problem subjected to various loading conditions (body force and 
thermal loading) was discussed Chauhan [28] has evaluated the 

5.1 F. using Transmission Photo- elasticity for the case of 
internal pressure loading He has confirmed that as the number 
of cracks increases, the S.I F. decreases. 

For the case of annular disc with large number of 
radial cracks subjected to centrifugal loading, no corresponding 
experimental results are reported in the literature. In this 
chapter, the preliminary studies conducted to evaluate S.I F. 
using Strain gauges and Photoelastic coatings are discussed. 

6.2 Determination of S.I.F throu^ Strain gauges : 

6.2.1 Strain field at the tip of a crack: 

It is shown in Ref [30] that a modified form of the 
Westergaard equations can be applied to describe the stress field 
associated with a single-ended crack in terms of two stress 
functions Z & Y by 

Oxx = KeZ - yImZ’ - yImY’ + 2ReY 
cTvv = + yImZ’ + ylmY’ 

UTxy = -yReZ’ - yReY’ - ImY 


. .( 6 . 1 ) 



where 


N t, 

Z (z) = S 2"-*^ 

n=0 
M 

Y (z) - 2 Bni2 , 2 = X + iy, 2 ’ = dZ/dz, Y’ = dY/dz 
m=0 

Substituting Eq (6.1) into the stress-strain relations gives, 

E ^xx = (1 ~ v)ReZ - (1 + v) yImZ’ 

- (1 4 V) yImY’ 4 2 KeY 

E €yy = (1 - v)ReZ 4 (1 4 V) ylmZ’ 

4 (1 4 V) yImY’ - 2vReY 


W-^xy = " yReZ’ ~ yReY’ - ImY ...(6 2) 

Of particular importance is the radial strain 
because this strain can be measured effectively at a large number 
of points in the field with commercially available strip gauges 
with 10 sensing elements on each strip. The radial strain is 
obtained from Eq.(6.2) by using the strain equation of 
transformation 

^rr = €xxCos*e 4 €yysin*e 4 Xj^ySinS cosS ...(6.3) 
which can be reduced to the following by using Eq. (6.2) and its 
derivatives. 


8 1 38 

2)i€j.y - Ac,n (kcos - sin8 sin cos2e 

2 2 2 

39 

4 sin® 6 cos8 cos ) 4 Bo{k 4 cos29) 

2 

3.8 9 

4 Ai r^ cos (k 4 sin* cos26 - sin® 8 cos9) 

2 2 

3 /2 

4 Bj r COS0 (k 4 1 - 6sin*9) 4 A 2 (k cos ) 

2 

3 e 

- sin cos2e - 3 sin®0 cos9 cos* 9) 

2 2 


4 2 B2(r*/(l4v))(l - (3 4 2v) sin*9 cos*e 


4 v sin* 8 (1 4 sin* 9)) 
where k = (l-v)/(l4v) 


. . .(6.4) 
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6.2.2 Least. Square Method to Evaluate S.I.F. ; 

In order to evaluate the strain field one can use strip 
gauges and these could be pasted along various S’s 

For ’n’ strain gauge sensors, Eq (6 4) may be used to 
form a system of equations in the unknown coefficients A^, Ai,A2) 

go,Bi and B 2 

o 1 

= Ao ri 2oi Bo ri hoi + Ai ri zn + Bi ri hn 
3/2 2 

+ A2 ri Z21 + ^2 ^21 

o _ ^ ^ ° ^ 1 3/2 2 

2w€rg - AornZon+Bornhon+Airn2in+Birnhin+A2rn22n+B2rnh2n 

.. (6 5) 

where g and h are functions of k and 8 given in Eq (6 4) 

Writing Eq. (6 5) in matrix form gives, 

{c} = [Dl {AB> ... (6 6) 

Since the total number of strain-gage readings exceed the 
number of unknown coefficients , i e., n > 6, Eq (6.6) represents 
an over-deterministic system of linear equations in the 
unknown coefficients. It is therefore necessary to statistically 
average the solution to make use of all the relevant data 
For this reason Eq (6 6) will be solved in the least-square 
sense For the least-squares solution of Eq (6.6), it is 
necessary to find a set coefficients {AB} which minimizes the 
vector {r> where 

{r} = {c} - [D3 {AB} .. (67) 

Ordinarily the solution to Eq (6 6) is obtained by forming the 
normal equations of the system, 

[D]'^ {c} = CB}'^ CD} CAB) 


. .( 6 . 8 ) 
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which yield the least squares solution for the unknown 
coefficients 


{AB} = ( [D]'^ [D3 [D]"^ {c} 


(6 9) 


Note that {AB} is unique for a given set of n data points When 
using a power series such as those given by Eq (6 l),the 
formation of normal equations frequently gives rise to 

numerical instabilities when the number of terms in the series 
becomes large Previous expeiience with boundary collocation 
[31] indicates that these instabilities occur when more than 15 
terms are used in the series expressions of Eq (6.1) 

Instabilities were not expected in this investigation since 
only SIX unknown coefficients were involved. 

A second method for the least-squares solution of Eq.(6.6) 
which avoids the instabilities associated with the formation of 
the normal equations is the OR decomposition. In this procedure 
a sequence of Householder transformations are applied to matrix 
resulting in the partitioned product 


[Q]D] 


m 

[o3 


. .( 6 . 10 ) 


where matrix [Q3 is orthogonal & matrix [R] is upper triangular 
If Eq (6.7) is multiplied by [Q] , 

£Q] {r} r: [Q]{c} - [Q3[D]{AB3 ...(6 11) 


and if the product [Q] {c} is partioned as [Q] {c} 


it may be rewritten as 


- 

{f} 

{g} 

. . . ( 6 . 12 ) 


{f} - CR3 {AB} 
{£} 


[Q3 {r} = 


. .(6 13) 
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The least-square solution to Eq (6 6) can then be found through 

standard elimination procedures on the determined system 

CR3 {AB} = {f} . (6 14) 

which minimizes the residual vector {r} Using the commercially 

available NAG subroutines, the least-squares solution to the 

linear system of Eq.(6.6) can be performed by OR decomposition 

A software has been developed to solve such set of 
equations using QR decomposition [32], 

6.2.3 Dse of a Single Strain Gauge to evaluate S.I.F. : 


In 

set to nil. 


Eq (6 4), contribution of Bq parameter 
if (k+cos2e) = 0 


can be 


( 1 - v) 

1 e cos2e = -k ..(6,15a) 

(1 + V) 

also, Ki = 48/3 XT ECj-j. . (6 16) 

Thus Ki can be determined by measuring the strain at 
a point close to the crack-tip on a radial line whose inclination 
is defined in terms of Poisson’s ratio of the material via 
Eq (6 15a), with another condition that the gauge must be 
placed at a position where radial co-ordinate r > 0.5 * t 

where t is the thickness of the specimen [33]. 


In case the orientation of the strain gauge (a) as 
shown in Fig. (6.1) is not the same as 6 , the inclination , then 
the governing equation for its placement becomes as 

tan 0/2 = -cot 2a ...(6.15b) 
For Aluminium, taking v = 0 3, 0 is found to be equal 

to 60“ . 

An experiment was conducted to verify this procedure to 
find stress intensity factor using one strain gauge using SEN 
specimen which is explained as below : 



6.2.4 Experiaental Details : 

Three specimens of size 310 * 35 * 3 2 mm were made of 

Aluminium and cracks were cut with a/w ratio equal to 01,02, 

0 3 respectively on the edges using 0.25 mm milling cutter. 

One strain gauge each was pasted at r - 3mm & 0 = 60® 

from the crack tip on both sides of the specimens of a/w equal to 

0 1 and 0 2 and at r = Burn & 6 = 60® lor the specimen of a/w 

equal to 0 3. Close-up view of one such strain gauge pasted is 
as given in Fig, (6 2) 

These strain gauges were connected in quarter bridge in 
the strain- indicator and the specimens were loaded in the MTS 
tensile testing machine with a grip length of 30 mm upto a 
maximum load of 5 kN in steps of 0 5 kN and the strains developed 
as recorded in the indicator were noted The experimental set-up 
IS shown in Fig. (6.3). 

The strain readings noted are as tabulated in 
Tab (8 1). Also tabulated are the percentage errors in S I F 
using strain gauge as compared to the S I.F. evaluated by 
boundary collocation method 

Although this method of strain-reading from single 
strain gauge pasted at a particular position and orientation 
requires the prior knowledge of the valid field and the plastic 
zone correction unlike the over-deterministic method, the results 
show that it could be be utilised to determine the Opening Mode 
Stress Intensity Factor (Kj) with accuracies sufficient for 
engineering applications 

Experiments were also planned to be conduci^^d on 
rotating discs to verify this method. But, as there was ylay in 
getting the slip rings required for the experiments, th# could 
not be conducted. 
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De-tenaination of S.I.F. through Birefringent-coatings : 

6.3.1 Introduction : 

Birefringent-coating method, a branch of photo- 
ticlty is a whole field techiiique to study the stress 
stribution on the surface of a body. In this method, a thin 
ayer of photoelastic plastic is bonded to the surface of the 
y under investigation. When load is applied to the body, the 
trains set up m it are transmitted to the photo elastic 
ing. The coating then becomes doubly refracting and a fringe 
pattern is visible when viewed in the field of polarised light 
birefringence observed is directly proportional to the 
intensity of strain developed The data obtained by this 
nique IS analogous to the data obtained from innumerable 
gauges of virtually zero gauge length A reflection type 
Po arii^cope is used to observe the birefringence. 

p. Consider a plane object with a thin coating as shown in 

^ ^ ^ E54] subjected to loading which produces a plane 

StSitO of V T • 

assumed that the stresses do not vary 
through the thickness of the coating. 

Let ajp, 02pj ^ip. and € 2 p refer to the principal 
^ resses and principal strains on the surface of the object and 
^lc« 02c> ^iQ & € 2 c correspond to the coating. If it is 
assumed that surface displacements are transmitted without 
distortion to the coating, then 

^3c = <^3p - Si 

!lc = €ip ..(6.17) 

^2g “ ^2p 

rom stress -optic law, the difference of principal 
stresses (cric - a2c) in the coating allowing for double passage 
of light is given by 

FoN 

- cr2c) = . .(6.18) 

2hc 

where N is the frixige order observed, 

F is the materials stress fringe value and 
h is the thickness of the coating. 
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Using Eqs.(6.17) , it can be shown 

Ep (1 + Yc) KFo 

~ <^2p'> = . (6.13) 

Ec (1 + Vp) 2hc 

where Ep & Vp are elastic modulus and Poisson’s ratio 
of the object and Eq & are elastic modulus and Poisson’s ratio 
of the coating 

Also, 

(^Ip ~ ^2p^ ~ NA/2h^K . . , (6.20) 

where K is the strain optic coefficient for the 
coating provided by the manufacturer and _X is the wave- 

length in nanometers of the light being used. 

Thus, from Eqs.(6.19), it is seen that the 

difference in principal stresses or principal strains can be 
determined from the isochromatic fringe order in the 
birefringent coating. Also, the directions of the principal 
stresses or strains are obtained by the isoclinios 

Isoclinic and Isochromatic are defined as follows • 

1. Isoclinic : The locus of points where the directions of the 
principal stresses coincide with a particular orientation of 
polarizer analyser combination . 

2. Isochromatic : It is the locus of points where the values of 

(o^ - 02 ) are such as to cause a relative phase difference of 2mx 

when the background is dark and (2m+l):ir , (m=0, 1, 2. . . ) when the 

background is bright 
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6-3.2 Expericnents conductied : 
i) S-E-H. Speclfflteut : 

A specimen of size 310 35 * 3 uua of Aluminium was 

prepared and a crack of length 7 mm (a/w=:0 2) was cut on one edge 
of the specimen using 0.25 mm milling cutter A photoelastic 
coating of size 50 >k 35 * 3 mm, of Type PS-IA (of Measurements 
roup Inc., U.o.A), of K factor .15 and thickness 3 mm was bonded 
onto the surface using PC-1 cement Note that a crack of the same 
size was also cut on the coating. 

The specimen was loaded in a MTS tensile testing 
machine gradually. The photograph of the fringes seen through 
the reflection polariscope at a load of 5,000 Newtons is shown as 
^ igs. (6. 5-6.6) . The fringes represent the stress field near the 
crack tip. The fringes show a mixed-mode loading (Mode-I + Mode- 
II) and this is due to the fact that crack cut was not exactly 
perpendicular to the loading direction 

ii) Sotating disc : 

A disc of size 260 OD, 30 ID, 3mm thickness was made of 
Aluminium Two internal cracks of length were cut diametrically 
opposite through 0.25 mm milling cutter and photoelastic coatings 
were bonded onto the crack tips. The disc was clamped and driven 
by a 1/8 HP AC motor The disc was rotated at different 
speeds through a Variac connected in series. 

The disc was viewed through reflection polariscope 
under stroboscopic illumination. The experimental set-up and 
the fringes seen at a speed of 3000 RPM are shown as 
Figs, (6. 7 6.8) However due to lack of time, this experiment 
could not be conducted in detail 
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6 . 4 Closure 

It is seen that both Photo-elastic coating and Strain- 
gauge technique could be employed to find SIP. Though Photo- 
elastic coating is a whole-field technique, the centrifugal 
stresses are not sufficient to produce sufficient fringes on the 
coating Use of a singe strain gauge to find S.I.F. is a 
promising technique and with modern strain gauge instrumentation, 
It is possible to find even 1 micro-strain precisely. Bence, with 
the use of slip rings, one can expect to get better estimates of 
S.I.F. by strain-gauge technique in rotating discs 





s. 

No 

Load 

(N) 

Micro 

Micro 

^Avg 

Micro 

HH 

Ml 

1 

500 

52 

68 

60 

65 

7 7 

2 

1000 

107 

134 

120 

130 

7 7 

3 

1500 

160 

204 

182 

195 

6 7 

4 

2000 

212 

270 

241 

259 

6 9 

5 

2500 

268 

342 

305 

324 

5 8 

6 

3000 

^ 320 

414 

367 

[ 

388 

5 4 

7 

3500 

* 376 

^ 482 

429 

453 

5 3 

8 

4000 

428 

1 

' 552 

490 

518 

5 4 

9 

4500 

482 

622 

552 

582 

5 2 

10 

5000 

540 

692 

616 

647 

4 8 


Table 8.1' Sbrain Gage Readings in SEN Specimen 
Strain Gage Length = 3 mm 
Strain Gage Resistance = 120 Q 
Strain Gage Factor - 2.0 
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FIG. 6, 4 Stresses in Photo-elastic Coating 





FIG. 6.6 Close-up View of Fringes seen in 
Birefringent-Coating 


CHAPTBK 7 


CC9ICL0SI0HS AND SDGGKSTKWS FOR FOTDRK WORK 


It has been shown that cyclic symmetry concept could be 
used with success to determine S.I.F in the case of a rotating 
annular disc with periodic radial cracks Considerable amount of 
work has been done to determine S.I.F’s as a function of crack 
lengths. number of cracks and b/a ratios for both, inner and 
outer cracks 

It has been further shown that the concept of cyclic 
symmetry can be used also to determine the S I F in the case of 
an annular disc subjected to Thermal loading Transient 
temperature distribution in the case of a clutch plate has been 
found out using a combination of FEM and FIW. 

Use of single strain gauge to find S.I.F is a 
promising technique with modern strain gauge instrumentation 
This has been verified finding S.I.F. of SEN specimen 

Following are suggested for future work ; 

1 The program developed could be modified 

incorporating Centrifugal force , Axial loading and Boundary 
conditions such as to simulate the actual conditions in a clutch 
plate. 


2 Strain-gauge technique could be utilised in 

finding S.I F in rotating discs using Slip rings 



MTmm h 


£HhU FUNClIONr. OF B-NODE ISOPARAMETRIC ELEMENT 

Ni = (1/4) n-r)(l-s)(-r-s-l) 

N2 = (1/4) (l+r)(l-s)(r-E-l) 

N 3 = (1/4) (l+r)(l+ 5 )(r+s-l) 

N 4 (1/4) ( 1-r) (1 + s) (-r+s-l ) 

N 5 = ( 1 / 2 ) (l-r^)(l-s) 

Ns = (1/2) (l-^r){l-s^) 

N 7 = (1/2) (l-r^)(l+s) 

Ns = ( 1 / 2 ) (l-r)(l-s^) 



APPENDIX - B 


tkmpkratohe distribution in a clutch plate of an automobile 

I> Introduction ■- 

Clutch plate of an automobile is subjected to 
centrifugal and thermal stresses when it engages with the engine 
flywheel to transmit torque to the transmission shaft. If the 
temperature distribution in clutch plate is known, corresponding 
stresses can be calculated easily using the software discussed 
Z J Jania [293 derives an analytical expression for the transient 
temperature distribution in such plates assuming constant rate of 
wear of the rubbing surfaces. With this assumption, the rate of 
energy dissipation and hence the temperature distribution is only 
a function of time and it is constant along the radii. But, if 
the pressure due to the axial clamping force is assumed to be 
constant over the rubbing surfaces, the distribution becomes a 
function of radii also An attempt has therefore been made in 
this Appendix to find the Transient temperature distribution 
which occurs in such clutch plates assuming constant pressure by 
combination of FEM and FDM teohnigues as described below ■ 

II) Numerical Formulation : 

Fig 7-1 shows two clutch plates in contact under the 
action of an axial clamping force F Plate II is coated with 
frictional material and plate I is made of steel. Surfaces AD & 
BC where rubbing occurs can be thought of as plane sources of 
heat producing kW per sq.mt. Temperature distribution in 

the Clutch plate I is governed by the Transient Fourier 
Conduction equation given below; 

15 dT Pc 6 T 

-- (r ) + 

r 6 r or bz 


k dt 


-..(7.1) 



where ’F’ ir Mass density, ’c’ is Heat capacitance and ’k’ is 
Thermal ct -nductivity of the material This equation is solved 
along with the Boundary and Initial conditions Before 
discussing the Boundary conditions, it is first necessary to 
estimate the Rate of Heat generation (qo(t)) 

i) Rate of Heat generation : 

A typical Power transmission can be idealised as shown 
in Fig 7 . 2 for the purpose of estimation of the rate of heat 
generation [29]. Here the Resistive load torque (Tr) is 
simulated by the brake For simplicity, ail torques are assumed 

constant 


When compliance of the system is neglected, 

Ti - Tc 

t + ri 


wi(t) = 


W2(t) 


Tc - T2 

T2 


t + r2 


(7.2) 


(7.3) 


where w(t) is angular velocity in rad/sec. , 

T is torque in N-m, 2 

I is Moment of Inertia in Newton-m-sec 
r is initial velocity in rad/sec , 
t xs time in sec. 

and subscripts 1 and 2 pertain to the 


^d load sides and ’c’ refers to the clutch 

At the end of slip period (to), 

wi(to) = V2itoy 
Therefore, t© ■" 


engine 


II l2 (ri - r2) 


Tc di t l 2 ) - (I2T1 + T1T2) 

It is assumed that the Pressure ’p’ due to the axial 
clamping force is uniformly distributed over the plate area. Then 
the Clutch frictional torque (dTc) over the elemental area 27trdr 

is given by, 


where u 


dTc = )ipr(2xrdr) 
is the coefficient of friction 


...(7 5) 


The rate at whach heat is generated over the element 
area 2Tirdr in the clutch during slipping is given by, 

QoCt) = jipr(wji (t)-W2 (t) ) . .(7 6) 

Thus, the rate of heat generation is directly 

proportional to the product of pressure and radial position 
Since the pressure is assumed to be constant, the rate of heat 
generation would be directly proportional to radial position of 
the plate. 

The effect of lubricating oil is not considered in 
calculation of the rate of heat generation as it is difficult to 
keep track of change in lubrication from hydrodynamic to boundary 
type while the oil is squeezed out of the space between the 
plates as they close Depending on the type of lubricant, under 
boundary lubrication conditions, temperature can be 30 to 60 
percent lower than in the case of dry friction. 

ii) Boundary conditions : 

At the rubbing interface AD & BC, the rate of heat 
generated is given by the Eqs (7 2-7 3 & 7. 5-7. 6). The ratios of 
heat flows into the plates I and II are not equal but depend on 
the thermal conductivity of the materials in contact If the 
thermal conductivities are k^ and Ik.2-, then the heat fluxes into 
the plates I and II given by , 

qi(t) = (k2/ki + k2) Qo 
q2(t) = (ki/ki + k2) q© 

Heat transfer from the clutch plates to the 
surroundings during slipping is assumed to be neglible 
Therefore, at the boundaries AB and CD, 

q(t) = 0. 

For axisymmetric problems only half the plate needs to 
be considered. Therefore the conditions are to be specified on 
the boundary IF also Here also it is assumed that q(t) = 0. 

Then at boundaries BB, 

q(t) - 0 .(7 8) 

Eqs (77 & 7.8) are the boundary conditions necessary 

for the solution of Eq.(7.1). 
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iii) Initial Conditions : 

In the beginning of engagement, the clutch temperature 
is the same as the ambient temperature i e , 

At time ’t’ = 0, T = T 0 ,(79) 

III) Solution Technique : 

i) GalerMn Integral : 

The first step in the Numerical solution of the 
Eq (7 1) is the formulation of the Galerkin integral This 
integral is obtained by multiplying the Eq.(7.1) with the weight 
function ’W’ and integrating the product by parts over the whole 
domain. Thus the expression for the integral becomes, 


I(Ti W) 


oT 6T oW oT 6W 

(pc V? + k (-- : — + ) 

5-t or or or oz 


27 irdrdz 


- j q W 2%r ds . . . / IV3 

where q* is the prescribed heat flux on the boundary T . 
Here, W is required to satisfy homogeneous version of the 
essential boundary conditions. 


ii) Finite Eleaent Approximation : 


In Galerkin Finite Element Method, the Temperature and 
the Weight function are approximated using same shape functions 


Thus over a typical element, 

T = {N}"^ {T}® 

W = {N}"^ {W}® 

where {N> is the Shape function vector and {T}® and {W}® 
the nodal values of T and W. 


. . (7 11) 
contain 



Zft 


Substituting Eq 7.11 into Eq 7.10, we get 

^e aT e e e e b (1 12) 

I = 2 {W}^ {C}^ {T}® + m CK3® {T}® - 2 {W}^ {Q}^^ 

e=l b=l 

where = number of area elements. 

n^ = number of line elements on the boundary T 

on which ’q’ is specified and 

pc[N]^ [N] 27rrdrd2 . 7.13 

= Elemental Heat capacity matrix 
T 

k EB3® [B3 2 xrdrd 2 . . .7.14 

= Elemental Heat conductivity matrix 

■> 

[Q3^ = Q* tN3^ 27[rdrdb 7 -lb 

i 




= Elemental Heat flux vector 
Writing I in terms of (liobal Matrices, 

I = {W}"^ [C3 {T} + {W}”^ [K3 {T} - {Q} . .7.16 

where [C3 , [K3 & {Q} are obtained by assembling the appropriate 
elemental matrices Since {W} f {0}, 

[C3 {T} + [K3’^ {T} - {Q} . .7.1' 

The time integration is done by Finite Difference 
Method as explained below 

iiDTinte Integration by Finite Difference Method : 

The time integration of the Eq.{7.17) is done using a 
Finite Difference scheme We approximate the Time derivative {T} 
by the following difference formula; 


e {T}i+i + (1-e) {T}i 


{Tli+i - {T>i 

(^t)i+l 


2 

+ 0(At)i+i 


where 0 takes the values as 0,1/2, 2/3 or 1 depending on 


.(7.18) 

whether 
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it IS I‘''iv}fird Crack~Nicolson, Galerkin, ot Backward 

Difffjence sch<=^5D^,* used Here, Backward difference scheme was used 
(0 = 1 ) 


Substituting Eq.7.18 into Eq.7 17, 

CA]a+i {T}ii = ^ 1 = 0^1^ 


wliere, 


(7 19) 


LAlx+1 - [c] + (^t)i+i e [k] . .(7.20) 

= {[c] - (^t)i+i (1 - 9) [k]} {T}i 

+ (>it)3^41 { 8{Q}i+l + (1 - 0) {Q}i } 

. (7.21) 

This equation is solved from the first time step upto 
the required time. Time limit was chosen as ’t’ = ’to> 

A sample problem as described below has been solved 
using the above mentioned method and the results are shown in 
Figs ( 5 1-5 3) . 


IV) Sample Problem : 

Consider a system in which the clutch is the multiple 
disc type with alternate steel and cork-coated plates System ajid 
clutch data are as below : 


II 


1 729 

Newton-m-sec 

l2 

r 

5.186 

Newton-m-sec 


Z 

346 

Newton-m 

T2 

Z 

276 

Newton-m 


z 

500 

Bad. /sec 


= 

200 

Rad /sec 

21 

z 

0.002 

m 

To 

i: 

00 

0 



z 

23.1 

W/m^-^F/m 

^2 

r 

0.231 

W/m^-“F/m 

Pi 


76,800 

3 

Newton/m 

Cl 


32.56 

Joule/Newton- 

refers 

to 

the clutch steel ■ 


where subscript ’1’ 

’2’ refers to the cork-coated friction plate 


plate 


S8 


and 



<7 V 



Friction Lining 


FIG, 7.1 Rubbing surfaces of a clutch plate 
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